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Abstract 



> . 

^ ■ In order to study light unstable nuclei systematically, we propose a new method 

. "AMD + Hartree-Fock" . This method introduces the concept of the single particle 

\ orbits into the usual AMD. Applying AMD + HF to Be isotopes, it is found that 

\ the calculated lowest intrinsic states with plus and minus parities have rather good 

On ' correspondence with the explanation by the two-center shell model. In addition, 
by active use of the single particle orbits extracted from AMD wave function, we 
construct the first excited 0^ state of ^'^Be. The obtained state appears in the 

O \ vicinity of the lowest 1~ state. This result is consistent with the experimental 

^ \ data. 

> 

X 
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I. INTRODUCTION 



For the study of the structure of hght unstable nuclei which are and have been extensively 
studied experimentally by the use of radioactive beams [1-3], various theories [4-20] have been 
used [21,22] . Among those theories, the antisymmetrized molecular dynamics (AMD) approach 
has already proved to be useful and successful [^ . One of the large merits of the AMD approach 
is that it does not rely on any model assumptions such as axial symmetry of the deformation, 
existence of the clustering and so on. 

Until now this AMD approach has not explicitly utilized the concept of the single particle 
motion in the mean field. The experimental data, however, show that very often we can get 
better understanding of the structure of unstable nuclei in terms of the single particle motion in 
the mean field. In order to see this point in more detail, we here discuss briefly some features 
of Be isotopes as an example. In Be isotopes there is the famous problem that the ground state 
^^Be has an anomalous parity. According to the usual shell model, its spin-parity should be 
I . But in experiments it is An answer to this problem is that due to the deformation the 
fi'^=| level in sd shell comes down below the upper ^'"=2 and the last neutron occupies the 
lowered il'^=^~^ . ^^Be is also famous as a halo nucleus. Since the lowered orbit contains s-orbit 
component, we can understand that the halo property is due to the long tail of the s-orbit. For 
the convenience of later discussion we call this lowered ^1'^=^^ level " halo level". 
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FIG. 1. The complete neutron level scheme of the two-center shell model as a function of the 
eccentricity parameter Zq. On the ordinate at the left are given the spherical shell model states of the 
original nucleus together with the energy scale, while on the right-hand side the quantum numbers of 
the spherical shell model states in the fragment nuclei are seen. This figure is quoted from [24]. 
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In order to discuss other Be isotopes than ^^Be, we use the single particle diagram by the 
two-center shell model [18-20], which is shown in FIG.l [^]. According to the AMD calculation 
for Be isotopes, we can support the idea that Be isotopes have the core part with an approximate 
dumbbell structure of two a clusters, which is the reason why we use the two-center shell model. 
^Be is known to have three rotational bands which have the band head states with J'^=| , | 
and We can explain these three bands in the two-center shell model as follows. We put the 
four neutrons in order from the lowest level at about Zo=2.5 fm in this model, and when the last 
neutron is put in the levels of f^''=| (a level), | (/3 level) and |^(7 level) , J'^ becomes | , | 
and I"*", respectively. Here, Zq is the eccentricity parameter, namely the distance between two 
centers in FIG.l. ^''Be is known to have two rotational bands which have the band head states 
with J'^=0"'" and 1~. In the same way as above, the four neutrons are put in order from the lowest 
level at about Zq=3 fm and when the configurations of the last two neutrons are [Q"^ = | ]^ and 
[^'^ = I Y[fl'^ = 1^]^, J'^ becomes 0"*" and 1~, respectively. Here we note that the Q'"= i"*" level 
is the halo level which appeared in ^^Be. 

These arguments show that by using the single particle orbits and by constructing the ex- 
cited state by particle-hole excitation we can study the structure of Be isotopes systematically. 
Therefore we recognize the importance of the single particle motion in the mean field. 

In this paper we propose the "AMD -|- Hartree-Fock" method (AMD-I-HF) which introduces 
the concept of the single particle orbits into the usual AMD. We extract the component of the 
mean field from the AMD wave function. Once we get the single particle wave functions in the 
mean field, we can study the structure of the excited states by constructing various particle-hole 
configurations. When we need elaborate description of the single particle motion like in the case 
of neutron halo phenomena, we need to improve the AMD single particle wave function in such 
a way that we adopt the superposition of several Gaussians instead of the single Gaussian wave 
packet in the usual AMD. 

This paper is composed of six sections. In the second section we explain the formalism of 
AMD-I-HF method. In the third section we show the single particle levels of the ground states of 
Be isotopes which are calculated with AMD-I-HF. In the forth section we show the single particle 
levels obtained with AMD-I-HF as a function of the deformation parameter. In the fifth section, 
based on the motivation of the AMD-I-HF, which means active use of the single particle levels, 
we construct the first plus-parity excited state of ^°Be by putting two neutrons into the halo level 
and study it. And the last section is for the discussion and the summary. 

II. AMD-HHF METHOD 
A. Variation under the constraint of deformation 

1. Hamiltonian and wave function 

In this paper, we used the Hamiltonian and wave function as described below. 
The Hamiltonian has the form: 

H = f + V, + VLs-fG. 
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Here T, Vc, Vls and Tq stand for the kinetic energy, the central force, the LS force and the 
center-of-mass kinetic energy, respectively. We neglected the Coulomb force, we used the Volkov 
No.l force p3 as the central force and the GSRS force PB| for the LS force. These forces have 
the following forms. 



where m=0.56, w=l-m, b=h=0, Vij=144.86 MeV, Va=-83.34 MeV, tr=0.82 fm, ta=1-Q fm, 
Vl5=900 MeV, ri5/j=0.20 fm, ri5^=0.36 fm. P(^0) is the projection operator onto the triplet 
odd state. 

The wave function of total system is expressed by a Slater determinant: 

1$ > = det[|(^i(j) >], 
< r\cp, > = (^)3/4^C7f exp[-Kr - |A >, 



where \Pi > is the spin-isospin function. In AMD+HF, single particle wave function is in general a 
linear combination of several Gaussians in order to describe single particle state more adequately 
than usual AMD. However in this paper we represented the single particle state by a single 
Gaussian wave packet for simplicity. |$ > is projected to parity eigenstates > and the 
energy variation is made after parity projection. 

1$=*= > = ^[|$ > ±P|$ >]. 
v2 



2. Cooling equation with constraint 

As mentioned above, the wave function is parameterized by Z"}. These parameters are 
determined by solving the following cooling equation, 

_ < ^^\H\^^ > 
~ < $±|<l>± > ■ 

Here {Xi} means {Cf, Z"}, A is an arbitrary real number and /i is a negative arbitrary number. 
We can easily prove that the energy of the total system decreases with time, 

dt 

We often need to obtain the minimum-energy state under some condition. For example, the 
condition is that the center-of-mass is fixed to the coordinate origin, or that the deformation 
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parameter is fixed to some value, etc. In such case, the condition is combined to the coohng 
equation by the analogy of the Lagrange-multiplier method. When the condition is represented 
as W(X°,X^*)=0, the previous coohng equation is changed as below: 

where rj is a, Lagrange-multiplier function, which is determined by = 0: 

dt 

^.dw,,, . .1 ,(M aw,-, 

From this equation we get 



•^ = E[(i-H77^7^ + (i + ^- 



^I'dX.dX* ' n'dX*dX,'' 

aw aw 



ax, ax; 

One may think that this method is applied only to the case where the number of conditions is 
one. But we can use it also in the case where several conditions exist. When several conditions 
are represented as 

Wi = 0,W2 = 0,---,W„ = 0, 

these conditions can be represented by a single equation, 

W = ci|Wi|' + --- + c„|W„|' = 0. 

Here ci ■ ■ ■ Cn are positive coefficients which adjust the difference of scale between Wi, W2, ■ ■ ■ , Wn. 
By this way, we can always use Eq.(l) without introducing several Lagrange- multipliers. 

3. Deformation parameter 
The deformation parameter (/?,7) can be obtained from the following relations. 
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where 



/ 2 ^ <^\iJ:iX^\^> 

< X > ^ ^^.^ , etc. 

< $1$ > 

Here x, y, z directions are the directions of the principal axes of inertia which are the coordinate 
axes of the body-fixed frame. We can easily confirm that this definition of deformation parameter 
is the same as that of Bohr model to the first order. Here expectation value of every operator is 
calculated with the intrinsic wave function before parity projection. In addition < >, <y^ > 
and < z"^ > is redefined so that < z"^ >>< y'^ >>< x"^ >. We must note the fact that in this 
definition z axis is not always the axial symmetric axis of the system. In the case of prolate 
deformation, z axis determined by this definition is certainly the axial symmetric axis, while in 
the case of oblate deformation it is not the axial symmetric axis. 

When we execute the cooling calculation with deformation constraint, it consumes time that in 
every step of the cooling we calculate the inertial axes to determine the deformation parameters. 
Therefore we use rotational invariant scalar quantities for calculating the deformation parameters 
so as to avoid the time-consuming procedure. We define the tensor quantity Q in any coordinate 
system as below: 

Qij — ^i^j ^1 

where i, j are x, y, z. Then calculating the trace of the tensors Q and QQ, we get the rotational 
invariant quantities TrQ and TtQQ. Since calculating these quantities in the space- fixed frame 
is identical to that in body-fixed frame, we obtain the following relations: 

TtQ space— fixed — TvQfjoijy— fixed 

= I] < >body-fixed= [< >< >< >] ^ ^ 
i i 

= [< x^Xy^Xz^ >]i. (3 + + ...), 

47r 

TrQQ space- fixed TrQQbody —fixed 

= E < >lody-fixed= [<X^><y^><Z^>]'^Y.'^t 
i i 

= [< X / X ^2 >] i • (3 + + . . .). 

TT 

Using these quantities, we define D as below: 
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= in space — fixed frame. 

Then we can calculate the deformation parameter /3 as 

Thus we can obtain /? without using the inertial axes. This approximation is enough correct to 
the second order of /3, which in ordinary case is less than unity. 

When wc impose the constraint that /? is fixed to some /^o, wc use D and Do corresponding 
to P and /So in the cooling calculation. The constraint condition function Wd is given as 



B. Single particle orbits described by AMD 

1. Method of calculating single particle orbits 

As mentioned in the introduction, we need to calculate the single particle orbits. We calculate 
them by the method shown below. This method extracts the single particle orbits from AMD 
wave function by mimicking the Hartrcc-Fock theory. 

In the Hartree-Fock method, the matrix elements hij of the single particle Hamiltonian in 
orthonormal base {\(f)i >} are given as 

hij = < 4>i\i\4>j > +X) < 4>i<t>k\v{\4>j4>k > -\<Pk(f>j >)• 

k 

By diagonalizing this h^j, we obtain the set of the single particle wave functions and the single 
particle energies. Now we construct Hartree-Fock-like single particle orbits and levels from total 
wave function of AMD. The single particle wave functions {cpi} of AMD are not orthonormal. 
So first we construct the orthonormal base by making linear-combination of {(fii}. We calculate 
the set of eigenvalues and eigenvectors of the overlap matrix B, 

3 

Bij = < (pi\(pj > . 

Eigenvectors {cf } are normalized. We define new base {|/" >} as follows: 

I/a > ^^J2cf\ipi > . 

We can confirm easily that this set forms an orthonormal base. 

Using this base |/" > we construct the matrix {ha/s} corresponding to the single particle 
Hamiltonian of Hartree-Fock as follows: 
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Kp = < fa\i\f(3 > + E < fafjlHlff^fj > -\flU >)■ 

7 

Diagonalizing this matrix {/la/j}, we calculate the eigenvalues {e^}, and eigenvectors {g^}-. 

The eigenvectors {g'^ are normahzed 
single particle energy is given as 

\P> 



We call such single particle orbits {\p >} and energies {e^} AMD-HF orbits and energies, respec- 
tively since they arc extracted from AMD wave function. 

Here we should note that these orbits and energies of AMD-HF are not totally equivalent to 
those of Hartree-Fock. The reason is as below. In AMD-I-HF, the Hartree-Fock equation is solved 
only within the functional space of single particle wave functions which is spanned by {i^i}. Thus 
in AMD-I-HF the Hartree-Fock self-consistency is satisfied only within this restricted functional 
space. We can regard the AMD-I-HF method as being a kind of restricted Hartree-Fock method. 
However this restriction has a large merit because the functional space spanned by the single 
particle wave functions {(^i} is obtained by the energy variation including parity projection. 



. The single particle wave function \p > belonging to the 

a 

a V i 

^. \/ LL 



2. Identification of AMD-HF orbits 

We need to identify the AMD-HF orbits \p > defined in the previous section. As almost 
nuclei seem to have axially symmetric deformation approximately, we expect that z component 
of the total angular momentum, j^, is approximately a good-quantum number, but that the 
magnitude of it, j^, is not. In addition, when the states that have the same absolute value of 
jz are degenerate, the expectation value of doesn't give us useful information. Therefore we 
use the square of ,31, so as to avoid this difficulty. For example, let us consider the state \p > 
which includes \ jz — m> and \ jz — —m > as below: 

I 1 r, , 1 

\p> ~ I [ \'m > +c\ — m >\. 



c 



2 



The expectation value of jz is 



IcP 



< P\jz\p > =m - - 

1 + c 
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On the other hand, the expectation value of jl is 

<p\3l\P > = 

from which we can identify m certainly. 

By the way, this identification is done on the body-fixed frame as noticed in the section 
about calculating the deformation. Here we calculate the inertial axes explicitly by diagonalizing 
the tensor Qjj, which was introduced in section II. A. 3. We choose the z-axis of the body-fixed 
frame to the direction of the Q^^^-'^^s eigenvector which has the largest eigenvalue. Here Q^p"^ 
represents Q calculated in the space-fixed frame. Using the eigenvector c*, we can relate 'fv"^^ 
and as follows: 

"Ahody _ q \space 
Jq — ^ J ) 

where j^P'^^'^ and the total angular momentum operator in space-fixed frame and that 

in the body-fixed frame, respectively. Thus the square in the body-fixed frame is represented 
using the quantities in the space-fixed frame as below: 

^ P\jz\P '^body— fixed ^ ^i ^ P\jijj\P space— fixed ■ 

ij 

In addition, we use also the square of the orbital angular momentum P, which is not a good- 
quantum number. This enables us to relate the obtained levels to the levels which would be 
obtained if the shape were spherical. 



III. AMD-HF ORBITS IN BE ISOTOPES 

In this section wc show the results about AMD-HF orbits of the ground states that are 
obtained by AMD. Here AMD single particle wave functions are composed of the single Gaussians. 

In TABLE 1, the AMD-HF levels are given. and indicate the expectation values of the 
operators and P respectively, which are obtained by the method mentioned in section 11. B. 2. 
(3 and B.E. are the deformation parameter and the binding energy respectively. We explain the 
results, about ^Be, ^°Be and ^^Be, that look interesting. Then we compare these with the simple 
model — two-center shell model — to help our understanding. 

a.^Be 

In the minus-parity AMD state of ^Be, of the highest level of neutron is 1.260, which is not 
equal to any of 0.25(=(|)^), 2.25(=(|)^), etc. This means the deformation is not axial symmetric. 
But L^=2.995 implies that this level contains rather large amount of the Ops component. 

In the plus-parity state, the last neutron is in the level of f2 = |. In addition of this level 
is 4.884. As this value is rather large, we can suppose that this level consists largely of d-orbit. 
In the two center shell model given in FIG.l, the plus-parity level named 7 level in FIG.l comes 
down in energy when the distance between two alpha clusters is large. This is consistent with 
our result because the deformation parameter /? obtained in AMD for the plus-parity state is 
0.801, which is very large. 
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^Be(+) B.E.=-21.26 MeV /3=0.437 



s.p. energy 


1^ 


occ. 


plus 


proton 








-2.48 0.250 


2.077 


2 


7.4 


-24.47 0.250 


0.118 


2 


96.0 


neutron 








-29.51 0.250 


0.148 


2 


93.1 


7Be(-) B.E.=- 


30.48 McV 


/3=0.566 


s.p. energy D."' 




occ. 


plus 


proton 








-11.04 0.250 


2.083 


2 


5.1 


-28.88 0.250 


0.219 


2 


96.5 


neutron 








-11.15 0.250 


2.059 


1 


4.3 


-28.87 0.250 


0.224 


1 


96.5 


-32.04 0.250 


0.440 


1 


82.2 


8Be(+) B.E.=- 


-48.00 MeV 


/3=0.656 


s.p. energy f2^ 


1^ 


occ. 


plus 


proton 








-18.99 0.250 


2.189 


1 


2.3 


-19.58 0.250 


2.127 


1 


0.2 


-34.69 0.250 


0.345 


2 


100.0 


neutron 








-18.87 0.250 


2.203 


1 


2.7 


-19.55 0.250 


2.131 


1 


0.3 


-34.69 0.250 


0.345 


2 


99.7 


''Be(-) B.E.=- 


46.96 MeV 


/3=0.521 


s.p. energy 


1^ 


occ. 


plus 


proton 








-21.82 0.249 


2.098 


2 


0.7 


-39.03 0.250 


0.231 


2 


99.7 


neutron 








+1.56 1.260 


2.996 


1 


21.8 


-18.18 0.269 


2.147 


1 


2.8 


-21.82 0.249 


2.098 


1 


0.7 


-33.70 0.256 


0.297 


1 


99.1 


-39.03 0.249 


0.230 


1 


99.7 


8Be(+) B.E.=- 


-42.97 MeV 


^=0.801 


s.p. energy f2^ 


1^ 


occ. 


plus 


proton 








-22.48 0.250 


2.287 


2 


3.1 


-35.04 0.250 


0.494 


2 


99.6 


neutron 








+6.79 0.250 


4.885 


1 


79.8 


-18.07 0.250 


2.280 


1 


4.5 


-22.37 0.250 


2.320 


1 


4.1 


-31.45 0.250 


0.483 


1 


98.6 


-35.06 0.250 


0.475 


1 


99.4 



iOBe(+) B.E.= 


-53.11 MeV 


^=0.381 


s.p. energy Cf 




occ. 


plus 


proton 








-23.00 0.254 


2.049 


2 


0.5 


-43.30 0.254 


0.108 


2 


99.6 


neutron 








-7.47 1.241 


2.487 


2 


12.7 


-19.55 0.278 


2.162 


2 


3.2 


-37.72 0.257 0.100 


2 


99.7 


lOBe(-) B.E.=- 


-44.95 MeV 


^=0.637 


s.p. energy H"^ 




OCC. 


plus 


proton 








-24.65 0.250 


2.106 


2 


0.0 


-39.75 0.250 


0.304 


2 


100.0 


neutron 








-3.41 0.667 3.387 


2 


59.4 


-20.39 0.264 


2.328 


2 


1.4 


-35.23 0.254 


0.270 


2 


99.7 


"Be(+) B.E.= 


-47.25 


McV 


,3=0.505 


s.p. energy il- 


1^ 


oce. 


)Ius 


proton 








-25.97 0.250 


2.089 


2 


1.9 


-44.97 0.250 


0.140 


2 


99.4 


neutron 








+4.25 0.272 


4.463 


1 


97.1 


-6.74 1.273 


2.153 


1 


3.7 


-8.84 1.253 


2.543 


1 


12.3 


-17.75 0.250 


2.327 


1 


2.5 


-22.47 0.258 


2.264 


1 


6.2 


-35.43 0.252 


0.147 


1 


99.8 


-39.55 0.255 


0.134 


1 


99.8 


iiBe(-) B.E.=- 


-54.70 MeV 


/3=0.271 


s.p. energy 




occ. 


j 


proton 








-24.97 0.250 


2.034 


2 


1.2 


-48.62 0.250 


0.042 


2 


99.7 


neutron 








-2.52 1.314 


2.350 


1 


8.2 


-7.15 1.241 


2.315 


1 


7.5 


-10.75 1.388 


2.311 


1 


7.4 


-18.04 0.251 


2.144 


1 


4.7 


-21.55 0.259 


2.115 


1 


2.9 


-36.61 0.250 


0.051 


1 


99.7 


-42.60 0.253 


0.035 


1 


100.0 



TABLE I. Single particle levels of the ground state of each 
nucleus. "B.E." and are the binding energy and the defor- 
mation parameter, respectively, "occ." indicates the occupation 
number of the level. " plus" indicates the percentage of plus-parity 
component. We regard two levels as the same level whose energy 
difference is less than 100 keV. In such case, the difference in 
and 1^ is about 0.01. 
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b. iOBe 

The calculated plus-parity state of ^°Be has a good correspondence with the two-center shell 
model, if we think our highest level of neutron corresponds to the level named the a level in 
FIG.l which comes from iz=\ of Ops in spherical case. 

But on the other hand, the minus-parity state does not have a similarity with the two-center 
shell model. The most surprising point is that the only three spectra of neutron appeared in our 
calculation. In the two-center shell model every single particle level is the eigenstate of parity. 
So in the model four spectra should appear. The explanation of this result by AMD-I-HF is as 
below: " The single particle orbit (AMD-HF orbit) extracted from AMD wave function is the 
parity-mixing state. And the last two neutrons are degenerate in that state." Explaining this 
phenomenon in terms of the two-center shell model is as follows. At some distance between two 
a particles, the levels Q = | and Q = i, coming from Ops and Ods respectively, approach each 
other. As a result, these levels are mixed, and then we see the only one level in which two levels 
are mixed. Here we investigate the percentage of parity-mixing in the last AMD-HF orbit of 
neutron. This AMD-HF level, which is represented as |/ >, can be written as follows: 

1/ > = ei|+ > +e2|- >, 
e? + e^ = l, 

where |+ > and |— > indicate plus- and minus-parity components, respectively. Using the 
operator V which reverses the parity, we calculate the following quantity A: 

A=<mf>-e\-e\. 

Then we obtain the ratio as follows: 

^1 ~ 2 ' ~ 2 ■ 

Calculated results are e^=0.60 and £2=0.40. As we expected this level is a state with large 
parity- mixing. 

c. ^^Be 

In the calculated plus-parity state of ^^Be, the last neutron occupies what we call "halo level" . 
f2^=0.25 means and L^=4.463 of this level means that this level contains a large component 
of 0d5. We think that this result is similar to that of the two-center shell model with medium 

2 

distance between two a particles. The obtained value of the deformation parameter /3=0.504, 
which is medium, is consistent with the medium distance in the two-center shell model. 

However in the calculated minus-parity state of -"^^Be. wc don't find so good similarity between 
AMD-I-HF and the two-center shell model. As /3=0.271 is rather small, we try to compare the 
obtained result to that of the two-center shell model with small inter-o; distance. In that case, 
the last neutron should be put into the level named the (5 level in FIG.l which has and 
comes from Op i in spherical case. But according to AMD-I-HF calculation, of the last neutron 

level is 1.314, which means that O is not |. 

By the way, there is the famous problem in ^^Be, which is that the parity of the ground state 
is not minus but plus. In this calculation, we have not succeeded to reproduce this anomalous 
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parity just like the previous AMD calculation of Ref. pSf. We think that this failure comes from 
the insufficient description of the single particle orbits in the present calculation where we have 
assigned only one Gaussian for one nucleon. 



IV. AMD-HF ORBITS WITH CONSTRAINED DEFORMATION 

In this section, we show the AMD-HF orbits at various f3, which look like Nilsson diagram. 
The levels at each f3 are calculated by the method as mentioned in the section II. A. 2. As in 
the previous section, we show only the results about ^Be(— ), ^''Be(— ) and ^^Be(— ), which look 
interesting. We show the diagram of the neutron's AMD-HF levels of ^Be(— ), ^°Be(— ) and 
^^Be(— ) in FIG. 2(a), FIG. 3(a) and FIG. 4(a), respectively. As the behavior of the highest level in 
each diagram looks especially interesting, we show the variation of L^, Q"^ and the percentage of 
the plus-parity component of the highest level in FIG. 2(b), FIG. 3(b) and FIG. 4(b). In addition 
as the second and the third levels in ^^Be(— ) look interesting too, we show the same quantities 
about these levels in FIG. 4(c) and FIG. 4(d). Notice that the scale on the right hand side is used 
for the percentage of the plus-parity component. 

a.9Be(-) 

n100% 

-75% 
-50% 
-25% 



0.2 0.4 0.6 0,8 1 02 0.4 0.6 0.8 T 

P P 
(a) diagram of AMD-HF levels (b) quantum numbers of the highest level 

FIG. 2. (a) Diagram of the AMD-HF levels of ^Be(— ). (b) Variation of L^, and the percentage 
of the plus-parity component of the highest level. 

In FIG. 2(a) we can see the inversion or the crossing of the single particle levels like that 
appearing in the two-center shell model. When we trace the behavior of the last neutron's level, 
in FIG. 2(b) we see that up to (3=0.67 the last level seems to be = | coming from Ops 
in spherical case. But beyond (3=0.67, its behavior changes apparently. Suddenly becomes 
rather large, becomes about 0.25 and the percentage of the plus-parity component becomes 
about 100%. Therefore we can recognize it as the lowered halo level that is = | coming from 
Ods. We notice that the behaviors of the other levels than the highest level are also changed 
beyond (3=0.67 in FIG. 2(a), but that L'^^Vl? and the percentage of the plus-parity component are 
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not changed. As the parity of the last neutron is almost plus beyond /3=0.67, the total intrinsic 
wave function is an almost plus-parity state. Therefore the projection to the minus-parity state 
is made by picking up the very small component of minus parity. 

b.iOBe(-) 

In FIG. 3(a), we can see that the last neutron's level keeps to be parity-mixed even if /3 is 
changed. When we look at the behavior of the quantum numbers ft and L in FIG. 3(b), we find 
that as (3 becomes larger \Q\ approaches to | and L increases too. This result seems to indicate 
that this level has two components, one being 101 = ^ coming from Ods and the other being 

— I coming from Ops, and that as P becomes larger the component of \fl\ = | increases. 
In addition, according to the parity-mixing ratio, this insight seems to be correct, since as /3 
becomes larger the ratio of plus-parity component increases. 
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FIG. 3. (a) Diagram of the AMD-HF levels of i°Be(-). (b) Variation of L^, n'^ and the percentage 
of the plus-parity component of the highest level. 



c.iiBe(-) 

We can see the tendency that the halo level, which is 101 = ^ coming from Ods, is lowered 
at large p. Especially at p—0.7 ~ 0.8 the tendency is prominent. In that region, we can see 

~ 2, O^ ~ 2 and the parity is 100% minus for the highest level in FIG.4(b). And in FIG.4(c) 
and 4(d), is rather large, O^ ~ 0.25 and the parity is about 90% plus for the second and 
third levels. Therefore we can identify in this region the highest level as |0| = | coming from 
Ops and the second and third levels as 101 = ^ coming from Ods which is the lowered halo level. 

2 I I Z 2 

Note that these two halo levels are not degenerate in this odd-neutron-number nucleus ^^Be. But 
beyond P—O.S the properties of the AMD-HF levels change again. Especially in FIG. 4(b) (the 
last neutron), the plus-parity component increases and becomes very large. We can explain 
this phenomenon with the two-center shell model as follows. The large deformation parameter 
P corresponds to the large inter-a distance in the two-center shell model. As shown in FIG.l, in 
such case the a level and 7 level approach very closely and are degenerate. In our calculation 
the same phenomenon as the two-center shell model seems to happen and then the component 
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of the halo orbit mixes into the highest level. 
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V. EXCITED POSITIVE PARITY BAND IN ^"BE 

The original AMD does not have the concept of the single particle orbit in the mean field. 
However we have confirmed in sections III and IV that AMD contains in it the mean field (and 
its single particle orbits). In sections III and IV, we used the single particle orbits only for 
"analyzing the AMD wave function". In this section, we truly "use" the single particle orbits by 
performing " particle- hole excitation" for the ground state. Our idea is that if the AMD-HF orbits 
are really meaningful, the description of excited states by the particle-hole excitation should be 
meaningful. 

Here we apply such p-h excitation for constructing the first excited 0+ state of ^°Be at 6.179 
MeV which is just above the lowest 1~ state at 5.960 MeV P^. One of possible configurations 



for this 0"^ is that two neutrons occupy the halo orbit, which is level coming from Ods in 

the two-center shell model. So we want to confirm that the second 0"*" constructed by such p-h 
excitation appears in the vicinity of the lowest 1~ as shown by the experiment. 

According to the two-center shell model, the lowest minus-parity state, which has the total 
angular momentum 1, seems to have the configuration that each of the levels, f2'^=f from Ops [a 

^ 2 

level) and from Ods (7 level) which is "halo level", is occupied by one neutron. Therefore 

by raising up the neutron occupying ^^^^=1 to VL^=]^ ^ the configuration becomes [fi'^ = 
which gives us the first excited plus-parity state having the total angular momentum 0. 

But according to the AMD-HF calculation of the lowest minus-parity state, we don't have 
good parity f^'^=| and f2'^=|^ levels but one parity-mixing level. In addition, the levels except 
this one are eigenstates of parity. Therefore in this case we can easily make the first excited 
plus-parity state from the lowest minus-parity state only by projecting out from the intrinsic 
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state of the lowest minus-parity state the plus-parity state: i.e. when the lowest minus-parity 
state |— g > is represented with its intrinsic state \^int > as 



\-9 > = > -T^l^int >], 

the first excited plus-parity state \+[ > is formed as 

l+i > = ;^[l*mt > +V\<^>int >]■ 

But this \+[ > may not be orthogonal to the lowest plus-parity state \+g > whose value of (3 is 
0.381. In order to avoid the mixing of this lowest plus-parity state, we orthogonalize > to 
> as below: 

+1 > = +; > — ' / • +9 > ■ 

Though 1+1 > and \±g > axe not the eigenstates of the total angular momentum, we can cal- 
culate the approximate energies of the projected states as below. For a given intrinsic state |$ >, 
the approximate expression for the energy Ej of the projected state with angular momentum / 
is 

^ r? 

Ei^< ^\H\^ > -—- < $1 J2|$ > +-=rl{l + 1), 
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In FIG. 5, we show the projected energies of ^°Be states as functions of (3. According to FIG. 5, 
the minimum energy and /? of each state is as below. 

lowest 0+ : -59.84 MeV at = 0.381 
lowest r : -51.90 MeV at /3 = 0.637 
second 0+ : -49.95 MeV at /3 = 0.852 

Therefore we can see that the first excited 0+ state appears in the vicinity of the lowest 1~ state 
with the small energy difference of 1.95 MeV. 

As mentioned previously, the AMD wave functions of Be isotopes have the core part with an 
approximate dumbbell structure of two alpha clusters. Therefore the inter-alpha distance is also 
useful as a measure of the deformation of Be isotopes. The above three values of /3 for 0^, 1~, 
and O2 in ^°Be correspond to the inter-alpha distance, 1.99 fm, 2.66 fm, and 3.55 fm, respectively. 

VI. SUMMARIZING DISCUSSION 

In this paper we have proposed a new theoretical method for the study of nuclear structure, 
which we have called the "AMD -|- Hartree-Fock" method (AMD+HF). The purpose of intro- 
ducing this new method is to develop the AMD approach by incorporating into it the concept of 
the single particle motion in the mean field. The new method consists of the following two steps. 
In the first step wc construct the AMD wave function by using usual AMD method. Then in 
the next step we extract the physical single particle motion in the mean field which is contained 
in the AMD wave function. This is done by diagonalizing the Hartree-Fock Hamiltonian in the 
functional space spanned by the single particle wave functions which constitute the AMD wave 
function. In this paper, for the sake of simplicity, we represented the AMD single particle wave 
function by the single Gaussian wave packet. But for the better description of the single particle 
motion we had better to represent the AMD single particle wave function by the superposition 
of Gaussian wave packets. Once we obtain the physical single particle wave functions in the 
mean field, we can utilize them for many purposes. First we can use them for understanding the 
physical situation of the ground state (or the ground rotational band states) of the system which 
the AMD wave function describes. Second we can use them for constructing the wave functions 
for the excited states of the system by constituting the particle-hole configurations. The latter 
way of utilizing the physical single particle wave functions in the mean field is beyond the scope 
of the ordinary AMD approach. 

The AMD-I-HF method is applicable for any nuclei, but it is expected to be especially useful 
for the study of unstable nuclei. The reason is as follows. The original AMD itself has proved 
to be very powerful for the study of unstable nuclei, which has been largely due to its basic 
character that it needs no model assumptions such as the axial symmetric deformation and the 
existence of the clusters. However, this merit of the original AMD means at the same time that 
we need to perform adequate analyses of the obtained AMD wave function in order to clarify the 
physics which the AMD wave function contains. As an example let us consider Be isotopes which 
we have treated in this paper. According to the AMD study of Be isotopes, the wave function 
of Be isotopes have the core part with an approximate dumbbell structure of two alpha clusters 
even for very neutron-rich isotopes near neutron dripline. Then there naturally arise questions 
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about what kind of dynamics governs so many neutrons distributed around the core part and 
about what kinds of interaction are existent between the core part and the neutrons. Our new 
approach, AMD+HF, can give us important information in answering to these questions. 

It should be noted that the AMD+HF method inherits most of the advantageous features 
of the original AMD. We need no model assumptions, we can perform variational calculation 
after projecting parity, we can use realistic effective nuclear force with finite range, we can 
perform angular momentum projection rather easily, we can superpose Slater determinants with 
no difficulty, and so on. Here we make a comment about the superposition of Slater determinants. 
In performing the ordinary AMD calculation, the superposition of Slater determinants causes no 



problem as was reported in previous papers [£3|. What we need to comment here is about the 



extraction of the single particle wave functions in the mean field contained in the superposed 
Slater determinants. In this paper, we have explained how to extract the mean field only in the 
case where the AMD wave function is expressed by a single Slater determinant. In the general 
case of the AMD wave function given by superposed Slater determinants, we first calculate the 
single particle density matrix from AMD wave function. Then we can calculate the Hartree- 
Fock type single particle Hamiltonian by the use of the density matrix. We will discuss this 
generalization elsewhere. 

As was pointed out in previous papers, the AMD method which adopts a single Gaussian wave 
packet for the single particle wave function is not suitable for the description of the long tails of 
the neutron halo orbits. For halo phenomena we have to adopt a superposition of Gaussian wave 
packets for the single particle wave function. However to represent all the single particle wave 
functions by the superposition of many Gaussian wave packets means a very heavy computational 
work. One of the aims of the AMD+HF method is to make more adequate description of the long 
tails of neutron halo orbits than the ordinary AMD method. By using the AMD+HF method 
we can identify the least bound neutron orbits which correspond to neutron halo orbits. What 
is necessary for us is to improve the single particle wave functions only for these least bound 
neutron orbits. Thus we expect that the AMD+HF method will enable us to treat neutron halo 
orbits in more efficient way than the usual AMD method. We will discuss this problem elsewhere. 

Generally speaking, the deformation of excited states is different from that of the ground state. 
Therefore in order to construct the wave functions with particle-hole excited configurations, 
we have to prepare the single particle wave functions for various magnitudes of quadrupole 
deformation. In this paper, we first calculated the minimum-energy AMD wave functions for 
various magnitudes of deformation by the use of the frictional cooling method under the constraint 
of the deformation and then extracted the single particle orbits from the obtained AMD wave 
functions. The investigation of the properties of the obtained single particle wave functions 
including the deformation- dependence of the single particle energies have shown that these AMD- 
HF orbits as functions of quadrupole deformation are rather similar to the single particle orbits of 
the two-center shell model quoted in section I. In spite of the overall similarity between AMD-HF 
orbits and two-center shell model orbits, there exist some interesting differences. Among them 
the appearance of the parity-mixed AMD-HF orbits is remarkable since it is out of the scope 
of the two-center shell model. A good example of the parity-mixed AMD-HF orbit is the least 
bound orbit extracted from the minus-parity AMD wave function of ^°Be. This orbit comes down 
in energy for large deformation and can be regarded as corresponding to the "halo orbit". 

In this paper we have studied Be isotopes by our new AMD+HF method. However the 
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application of the AMD+HF method to the study of excited states by constructing the wave 
functions with particle-hole excited configurations has been made only for one problem, namely 
the study of the second 0+ state of ^°Be at 6.18 MeV. According to the calculation reported in 
this paper, the O"*" state with the configuration of the two valence neutrons occupying the valence 
halo orbit which is parity mixed have the excitation energy near that of the lowest 1^ state. 
This result is consistent with the experimental data which show that the second 0"*" state at 6.18 
MeV exists in the vicinity of the lowest 1~ state at 5.96 MeV. We here would like to comment 
that this calculated result can be improved by using the density-dependent force. The density- 
dependent force works attractively at low density, and repulsively at high density. Therefore in 
the case where the density- dependent force is added to the force used in this paper, the second 
0"*" state comes down in energy more than the lowest 1~ state does, because the second 0"*" state 
is deformed more largely than the lowest 1~ state according to the calculated result. We expect 
that the second 0+ state appears more closely to the lowest 1~ state. 

Our study of Be isotopes in this paper was made for the region of ^Be ~ ^^Bc. But of course 
we can study ^^Be and ^^Be in the same way. Especially it is very interesting to study the role 
of the "halo orbit" not only in the ground state but also in the excited states in these ^^Be and 
i^Be. 

In summary, the AMD-HF theory can be expected to be not only useful for the analysis of 
the AMD wave function but also powerful for the systematic study of excited states in a different 
way from the ordinary AMD approach by relying on the concept of the particle-hole excitation. 
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